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1 Introduction 

Convex polyhedra, i.e., the intersections of finitely many closed affine half- 
spaces in R'', are important objects in various areas of mathematics and other 
disciplines. In particular, the compact ones among them (polytopes), which 
equivalently can be defined as the convex hulls of finitely many points in R'', 
have been studied since ancient times (e.g., the platonic solids). Polytopes 
appear as building blocks of more complicated structures, e.g., in (combina- 
torial) topology, numerical mathematics, or computer aided design. Even in 
physics polytopes are relevant (e.g., in crystallography or string theory). 

Probably the most important reason for the tremendous growth of interest 
in the theory of convex polyhedra in the second half of the 20"^ century was 
the fact that linear programming (i.e., optimizing a linear function over the 
solutions of a system of linear inequalities) became a widespread tool to solve 
practical problems in industry (and military). Dantzig's Simplex Algorithm, 
developed in the late 40's, showed that geometric and combinatorial know- 
ledge of polyhedra (as the domains of linear programming problems) is quite 
helpful for finding and analyzing solution procedures for linear programming 
problems. 

Since the interest in the theory of convex polyhedra to a large extent 
comes from algorithmic problems, it is not surprising that many algorithmic 
questions on polyhedra arose in the past. But also inherently, convex poly- 
hedra (in particular: polytopes) give rise to algorithmic questions, because 
they can be treated as finite objects by definition. This makes it possible 
to investigate (the smaller ones among) them by computer programs (like 
the polymake-system written by Gawrilow and Joswig, see ||2^ and j2^,^). 
Once chosen to exploit this possibility, one immediately finds oneself con- 
fronted with many algorithmic challenges. 

This paper contains descriptions of 35 algorithmic problems about poly- 
hedra. The goal is to collect for each problem the current knowledge about its 
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computational complexity. Consequently, our treatment is focused on theo- 
retical rather than on practical subjects. We would, however, like to mention 
that for many of the problems computer codes are available. 

Our choice of problems to be included is definitely influenced by personal 
interest. We have not spent particular efforts to demonstrate for each problem 
why we consider it to be relevant. It may well be that the reader finds other 
problems at least as interesting as the ones we discuss. We would be very in- 
terested to learn about such problems. The collection of problem descriptions 
presented in this paper is intended to be maintained as a (hopefully growing) 
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Almost all of the problems are questions about polytopes. In some cases 
the corresponding questions on general polyhedra are interesting as well. It 
can be tested in polynomial time whether a polyhedron specified by linear 
inequalities is bounded or not. This can be done by applying Gaussian elim- 
ination and solving one linear program. 

Roughly, the problems can be divided into two types: problems for which 
the input are "geometrical" data and problems for which the input is "com- 
binatorial" (see below). Actually, it turned out that it was rather convenient 
to group the problems we have selected into the five categories "Coordinate 
Descriptions" (Sect.||), "Combinatorial Structure" (Sect. ^), "Isomorphism" 
(Sect. 1^), "Optimization" (Sect. ||), and "Realizability" (Sect. ||). Since the 
boundary complex of a simplicial polytope is a simplicial complex, studying 
polytopes leads to questions that are concerned with more general (polyhe- 
dral) structures: simplicial complexes. Therefore, we have added a category 
"Beyond Polytopes" (Sect. |^), where a few problems concerned with general 
(abstract) simplicial complexes are collected that are closely related to similar 
problems on polytopes. We do not consider other related areas like oriented 
matroids. 

The problem descriptions proceed along the following scheme. First input 
and output are specified. Then a summary of the knowledge on the theoretical 
complexity is given, e.g., it is stated that the complexity is unknown ("Open") 
or that the problem is A/'P-hard. This is done for the case where the dimension 
(usually of the input polytope) is part of the input as well as for the case of 
fixed dimension; often the (knowledge on the) complexity status differs for 
the two versions. After that, comments on the problems are given together 
with references. For each problem we tried to report on the current state 
of knowledge according to the literature. Unless stated otherwise, all results 
mentioned without citations are either considered to be "folklore" or "easy 
to prove." At the end related problems in this paper are listed. 

For all notions in the theory of polytopes that we use without explanation 
we refer to Ziegler's book Similarly, for the concepts from the theory 
of computational complexity that play a role here we refer to Carey and 
Johnson's classical text . Whenever we talk about polynomial reductions 
this refers to polynomial time Turing-reductions. For some of the problems 
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the output can be exponentially large in the input. For these problems the 
interesting question is whether there is a polynomial total time algorithm, 
i.e., an algorithm whose running time can be bounded by a polynomial in the 
sizes of the input and the output (in contrast to a polynomial time algorithm 
whose running time would be bounded by a polynomial just in the input 
size) . Note that the notion of "polynomial total time" only makes sense with 
respect to problems which explicitly require the output to be non-redundant. 

A very fundamental result in the theory of convex polyhedra is due to 
Minkowski |Q and Weyl |6j] . For the special case of polytopes (to which we 
restrict our attention from now on) it can be formulated as follows. Every 
polytope P C R'' can be specified by an H- or by a V-description. Here, 
an Ti-description consists of a finite set of linear inequalities (defining closed 
affine half-spaces of M!^) such that P is the set of all simultaneous solutions to 
these inequalities. A V-description consists of a finite set of points in M!^ whose 
convex hull is P. If any of the two descriptions is rational, then the other one 
can be chosen to be rational as well. Furthermore, in this case the numbers in 
the second description can be chosen such that their coding lengths depend 
only polynomially on the coding lengths of the numbers in the first description 
(see, e.g., Schrijver Q). In our context, H- and V-descriptions are usually 
meant to be rational. By linear programming, each type of description can 
be made non-redundant in polynomial time (though it is unknown whether 
this is possible in strongly polynomial time, see Problem |2^). 

One of the basic properties of a polytope is its dimension. If the polytope 
is given by a V-description, then it can easily be determined by Gaussian 
elimination (which, carefully done, is a cubic algorithm; see, e.g., [^). If the 
polyhedron is specified by an H-description, computing its dimension can be 
done by linear programming (actually, this is polynomial time equivalent to 
linear programming). 

Furthermore, some of the problems may also be interesting in their po- 
lar formulations, i.e., with "the roles of H- and V-descriptions exchanged." 
Switching to the polar requires to have a relative interior point at hand, 
which is easy to obtain if a V-description is available, while it needs linear 
programming if only an "H-description is specified. 

We will especially be concerned with the combinatorial types of polytopes, 
i.e., with their face lattices (the sets of faces, ordered by inclusion). In particu- 
lar, some problems will deal with the k-skeleton of a polytope, which is the set 
of its faces of dimensions less than or equal to fc, or with its f -vector, i.e., the 
vector (/o(P), fi{P), ■ ■ ■ , fd{P)), where fi{P) is the number of i-dimensional 
faces (i-faces) of the d-dimensional polytope P (d-polytope) . Talking of the 
face lattice Cp oi a, polytope P will always refer to the lattice as an abstract 
object, i.e., to any lattice that is isomorphic to the face lattice. In particu- 
lar, the lattice does not contain any information on coordinates. Similarly, 
the vertex-facet incidences of P are given by any matrix (a„/) with entries 
from {0, 1}, whose rows and columns are indexed by the vertices and facets 



4 



Volker Kaibel and Marc E. Pfetsch 



of P, respectively, such that a^f = 1 if and only if vertex v is contained 
in facet /. Note that the vertex- facet incidences of a polytope completely 
determine its face lattice. 

A third important combinatorial structure associated with a polytope P is 
its (abstract) graph Qp, i.e., any graph that is isomorphic to the graph having 
the vertices of P as its nodes, where two of them are adjacent if and only if 
their convex hull is a (one-dimensional) face of P. For simple polytopes, the 
(abstract) graph determines the entire face lattice as well (see Problem p^ . 
However, for general polytopes this is not true. 

Throughout the paper, n refers to the number of vertices or points in 
the given V-description, respectively, depending on the context. Moreover, 
m refers to the number of facets or inequalities in the given 7i-description, 
respectively, and d refers to the dimension of the polytope or the ambient 
space, respectively. 

Acknov^rledgment: We thank the referee for many valuable comments and 
Giinter M. Ziegler for carefully reading the manuscript. 



2 Coordinate Descriptions 

In this section problems are collected whose input are geometrical data, i.e., 
the Ti- or V-description of a polytope. Some problems which are also given 
by geometrical data appear in Sections ^ and ^ 

1. Vertex Enumeration 

Input: Polytope P in W-description 
Output: Non-redundant V-description of P 

Status (general): Open; polynomial total time if P is simple or simplicial 
Status (fixed dim.): Polynomial time 

Let d = dim(P) and let m be the number of inequalities in the input. It is well 
known that the number of vertices n can be exponential (f2(mL'^/^J )) in the 
size of the input (e.g., Cartesian products of suitably chosen two-dimensional 
polytopes and prisms over them). 

Vertex Enumeration is strongly polynomially equivalent to Problem || 
(see Avis, Bremner, and Seidel j^]). Since Problem ^ is strongly polynomially 
equivalent to Problem || as well. Vertex Enumeration is also strongly 
polynomially equivalent to Problem ^j. 

For fixed d, Chazelle found an O^m^'''^'^^) polynomial time algorithm, 
which is optimal by the Upper Bound Theorem of McMuUen |Q. There 
exist algorithms which are faster than Chazelle's algorithm for small n, e.g., 
an 0(771 log n + (TOn)^^^/(L'i/2j-i-i) polylogm) algorithm of Chan 
For general d, the reverse search method of Avis and Fukuda ||^ solves the 
problem for simple polyhedra in polynomial total time, using working space 
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(without space for output) bounded polynomially in the input size. An al- 
gorithm of Bremner, Fukuda, and Marzetta |^ solves the problem for sim- 
plicial polytopes. Note that these algorithms need a vertex of P to start 
from. Provan js^ gives a polynomial total time algorithm for enumerating 
the vertices of polyhedra arising from networks. 

There are many more algorithms known for this problem - none of them is 
a polynomial total time algorithm for general polytopes. See the overview 
article of Seidel |Q . Most of these algorithms can be generalized to directly 
work for unbounded polyhedra, too. 

Related problems: |, |, |, 

2. Facet Enumeration 

Input: Polytope P in V-description with n points 
Output: Non-redundant 7i-description of P 

Status (general): Open; polynomial total time if P is simple or simplicial 
Status (fixed dim.): Polynomial time 

In iQ it is shown that Facet Enumeration is strongly polynomially equiv- 
alent to Problem || and thus to Problem |^ (see the comments there) . 
For this problem, one can assume to have an interior point (e.g., the ver- 
tex barycenter). Facet Enumeration is sometimes called the convex hull 
problem. 

Related problems: |l|, ^, ^ 

3. Polytope Verification 

Input: Polytope P given in 7i-description, polytope Q given in V-description 
Output: "Yes" ii P = Q, "No" otherwise 

Status (general): Open; polynomial time if P is simple or simplicial 
Status (fixed dim.): Polynomial time 

Polytope Verification is strongly polynomially equivalent to Problem |^ 
and Problem H (see the comments there). 

Polytope Verification is contained in coMV: we can prove Q ^ P by 
showing that some vertex of Q violates one of the inequalities describing P. 
li Q <Z P with Q ^ P then there exists a point p oi P\Q with "small" co- 
ordinates (e.g., some vertex of P not contained in Q) and a valid inequality 
for Q, which has "small" coefficients and is violated by p (e.g., an inequal- 
ity defining a facet of Q that separates p from Q). However, it is unknown 
whether Polytope Verification is in A/'P. 

Since it is easy to check whether Q C F, Polytope Verification is Prob- 
lem H restricted to the case that Q C P. 

Related problems: |l|, ^, ^ 
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4. PoLYTOPE Containment 

Input: Polytope P given in 7i-description, polytope Q given in V-description 
Output: "Yes" if P C Q, "No" otherwise 
Status (general): coAAP-complete 
Status (fixed dim.): Polynomial time 

Freund and Orlin proved that this problem is coA/'P-complete. Note 
that the reverse question whether Q C P is trivial. The questions where 
either both P and Q are given in 7i-description or both in V-description can 
be solved by linear programming (Problem |^) , see Eaves and Freund JlTf . 
For fixed dimension, one can enumerate all vertices of P in polynomial time 
(see Problem and compare the descriptions of P and Q (after removing 
redundant points). 

Related problems: ^ 

5. Face Lattice of Geometric Polytopes 

Input: Polytope P in H-description 
Output: Hasse-diagram of the face lattice of P 
Status (general): Polynomial total time 
Status (fixed dim.): Polynomial time 

See comments on Problem |l|. Many algorithms for the Vertex Enumer- 
ation Problem in fact compute the whole face lattice of the polytope. 
Swart ||60| , analyzing an algorithm of Chand and Kapur [^], proved that 
there exists a polynomial total time algorithm for this problem. For a faster 
algorithm see Seidel |5^. Fukuda, Licbling, and Margot |22| gave an algo- 
rithm which uses working space (without space for the output) bounded 
polynomially in the input size, but it has to solve many linear programs. 
For fixed dimension, the size of the output is polynomial in the size of the 
input; hence, a polynomial total time algorithm becomes a polynomial algo- 
rithm in this case. 

The problem of enumerating the fc-skeleton of P seems to be open, even if 
k is fixed. Note that, for fixed k, the latter problem can be solved by linear 
programming (Problem |2^) in polynomial time if the polytope is given in 
V-description rather than in 7i-description. 

Related problems: |l|, |, |, |l| 

6. Degeneracy Testing 

Input: Polytope P in Ti-description 

Output: "Yes" if P not simple (degenerate), "No" otherwise 
Status (general) : Strongly A/'P-complete 
Status (fixed dim.): Polynomial time 
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Independently proved to be A/''P-coniplete in the papers of Chandrasekaran, 
Kabadi, and Murty and Dyer 0. Fukuda, Liebling, and Margot ^ 
proved that the problem is strongly A/'T'-complete. For fixed dimension, one 
can enumerate all vertices in polynomial time (see Problem ^ and check 
whether they are simple or not. 

Bremner, Fukuda, and Marzetta [|| noted that if P is given in V-description 
the problem is polynomial time solvable: enumerate the edges (1-skeleton, see 
Problem ||) and apply the Lower Bound Theorem. 

Erickson ||l^ showed that in the worst case [2{m^'^^'^^^^ + to log m) sideness 
queries are required to test whether a polytope is simple. For odd d this 
matches the upper bound. A sideness query is a question of the following 
kind: given d + 1 points po, . . . ,pd in R'', does po lie "above", "below", or 
on the oriented hyperplane determined by pi,p2, • ■ • ,Pci- 

Related problems: |l|, ^ 

7. Number of Vertices 

Input: Polytope P in ^-description 
Output: Number of vertices of P 
Status (general): #'P-complete 
Status (fixed dim.): Polynomial time 

Dyer and Linial |^ independently proved that Number of Vertices 
is #7^-complctc. It follows that the problem of computing the /-vector of P 
is ^7^-hard. Furthermore, Dyer JT^ proved that the decision version ("Given 
a number k, does P have at least k vertices?") is strongly TVT^-hard and 
remains A/'T^-hard when restricted to simple polytopes. It is unknown whether 
the decision problem is in J\fV. 

If the dimension is fixed, one can enumerate all vertices in polynomial time 
(see Problem ^. 

Related problems: |l|, |lj 

8. Feasible Basis Extension 

Input: Polytope P given as {a; G R'' : Ax ^ b,x > 0}, a set 5 C {1, . . . , s} 
Output: "Yes" if there is a feasible basis with an index set containing S, 

"No" otherwise 
Status (general): A/'T^-complete 
Status (fixed dim.): Polynomial time 

See Murty |49| (Garey and Johnson Q, Problem [MP4]). For fixed dimen- 
sion, one can enumerate all bases in polynomial time. 

The problem can be reformulated as follows. Let P be defined by a finite set 
H of afRne halfspaces and let S" be a subset of H. Docs f]{H £ H : H ^ S} 
contain a vertex which is also a vertex of P? 
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9. Recognizing Integer Polyhedra 

Input: Polytope P in Ti-description 

Output: "Yes" if P has only integer vertices, "No" otherwise 
Status (general): Strongly coA/'T'-complete 
Status (fixed dim.): Polynomial time 

The hardness-proof is by Papadimitriou and Yannakakis | ]5l[ . For fixed di- 
mension, one can enumerate all vertices (Problem |^) and check whether they 
are integral in polynomial time. 

10. Diameter 

Input: Polytope P in ^-description 
Output: The diameter of P 
Status (general): AfV-hard 
Status (fixed dim.): Polynomial time 

Frieze and Teng |^ gave the proof of A/'T'-hardness. For fixed dimension, 
one can enumerate all vertices (Problem |l|), construct the graph and then 
compute the diameter in polynomial time. 

The complexity status is unknown for simple polytopes. For simplicial poly- 
topes the problem can be solved in polynomial time: Since simplicial poly- 
topes have at most as many vertices as facets, one can enumerate their ver- 
tices (see Problem 0), and finally compute the graph (and hence the diameter) 
from the vertex-facet incidences in polynomial time. 

If P is given in V-description, one can compute the graph (1-skeleton, see 
Problem ||) and hence the diameter in polynomial time. 



11. Minimum Triangulation 

Input: Polytope P in V-description, positive integer K 

Output: "Yes" if P has a triangulation of size K or less, "No" otherwise 

Status (general): TVT^-complete 

Status (fixed dim.): A/'T'-complete 

A triangulation T of a d-polytope P is a collection of d-simplices, whose 
union is P, their vertices are vertices of P, and any two simplices intersect 
in a common face (which might be empty). In particular, T is a (pure) d- 
dimensional geometric simplicial complex (see Section |7|) . The size of T is the 
number of its d-simplices. Every (convex) polytope admits a triangulation. 
Below, De Loera, and Richter-Gebert [^j^ proved that Minimum Triangu- 
lation is A/'T-'-complete for (fixed) > 3. Furthermore, it is TVP-hard to 
compute a triangulation of minimal size for (fixed) d > 3. 
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12. Volume 

Input: Polytope P in W-dcscription 
Output: The volume of P 
Status (general): #P-hard, FPRAS 
Status (fixed dim.): Polynomial time 

Dyer and Frieze [|l5j showed that the general problem is ^P-hard (and #P- 
easy as well). Dyer, Frieze, and Kannan | |T^ found a Fully Polynomial Ran- 
domized Approximation Scheme (FPRAS) for the problem, i.e., a family 
(^£)£>o of randomized algorithms, where, for each e > 0, computes a 
number with the property that the probability of (1 — e) vol(P) < < 
(1 + e) vol(P) is at least |, and the running times of the algorithms A^ are 
bounded by a polynomial in the input size and K 

For fixed dimension, one can first compute all vertices of P (see Problem ^ 
and its face lattice (see Problem ^ both in polynomial time. Then one can 
construct some triangulation (see Problem |l^) of P (e.g., its barycentric 
subdivision) in polynomial time and compute the volume of P as the sum of 
the volumes of the (maximal) simplices in the triangulation. 
The complexity status of the analogue problem with the polytope specified 
by a V-description is the same. 



3 Combinatorial Structure 

In this section we collect problems that are concerned with computing certain 
combinatorial information from compact descriptions of the combinatorial 
structure of a polytope. Such compact encodings might be the vertex-facet 
incidences, or, for simple polytopes, the abstract graphs. An example of such 
a problem is to compute the dimension of a polytope from its vertex-facet 
incidences. Initialize a set S by the vertex set of an arbitrary facet. For each 
facet F compute the intersection of S with the vertex set of F. Replace S by 
a maximal one among the proper intersections and continue. The dimension 
is the number of "rounds" performed until S becomes empty. 

All data is meant to be purely combinatorial. For all problems in this 
section it is unknown if the "integrity" of the input data can be checked, 
proved, or disproved in polynomial time. For instance, it is rather unlikely 
that one can efficiently prove or disprove that a lattice is the face lattice of 
some polytope (see Problems ^ . 

Sometimes, it might be worthwhile to exchange the roles of vertices and 
facets by duality of polytopes. Our choices of view points have mainly been 
guided by personal taste. 

Some orientations of the abstract graph Qp oi a. simple polytope P play 
important roles (although such orientations can also be considered for non- 
simple polytopes, they have not yet proven to be useful in the more general 
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context). An orientation is called a unique-sink orientation (US-orientation) 
if it induces a unique sink on every subgraph of Qp corresponding to a non- 
empty face of P. A US-orientation is called an abstract objective function ori- 
entation (AOF-orientation) if it is acyclic. General US-orientations of graphs 
of cubes have recently received some attention (Szabo and Welzl |^). AOF- 
orientations were used, e.g., by Kalai Since their linear extensions are 
precisely the shelling orders of the dual polytope, they have been considered 
much earlier. 



13. Face Lattice of Combinatorial Polytopes 

Input: Vertex- facet incidence matrix of a polytope P 
Output: Hasse-diagram of the face lattice of P 
Status (general): Polynomial total time 
Status (fixed dim.): Polynomial time 

Solvable in ©(minjm, n} ■ a ■ ip) time, where m is the number of facets, n is 
the number of vertices, a is the number of vertex- facet incidences, and ip is 
the size of the face lattice j33|. Note that (p is exponential in d (for fixed d 
it is polynomial in m and n). Without (asymptotically) increasing the run- 
ning time it is also possible to label each node in the Hasse diagram by the 
dimension and the vertex set of the corresponding face. 
It follows from ]33| that one can compute the Hasse-diagram of the k-skeleton 
(i.e., all faces of dimensions at most k) of P in 0{n ■ a ■ p-'') time, where (p-'' 
is the number of faces of dimensions at most k. Since the latter number is in 
0{n''+^), the fc-skeleton can be computed in polynomial time (in the input 
size) for fixed k. 

Related problems: ||, |lj 



14. /-Vector of Combinatorial Polytopes 

Input: Vertex- facet incidence matrix of a polytope P 

Output: /-vector of P 

Status (general): Open 

Status (fixed dim.): Polynomial time 

By the remarks on Problem it is clear that the first k entries of the 
/-vector can be computed in polynomial time for every fixed k. 
If the polytope is simplicial and a shelling (or a partition) of its boundary 
complex is available (see Problems f?] and |l^), then one can compute the 
entire /-vector in polynomial time [pSTChap. 8]. 

Related problems: H 
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15. Reconstruction of Simple Polytopes 

Input: The (abstract) graph Qp of a simple polytope P 

Output: The family of the subsets of nodes oiGp corresponding to the vertex 

sets of the facets of P 
Status (general): Open 
Status (fixed dim.): Open 

Blind and Mani proved that the entire combinatorial structure of a simple 
polytope is determined by its graph. This is false for general polytopes (of di- 
mension at least four) , which is the main reason why we restrict our attention 
to simple polytopes for the remaining problems in this section. Kalai ||35|| gave 
a short, elegant, and constructive proof of Blind and Mani's result. However, 
the algorithm that can be derived from it has a worst-case running time that 
is exponential in the number of vertices of the polytope. 

In [^2| it is shown that the problem can be formulated as a combinatorial 
optimization problem for which the problem to find an AOF-orientation of Qp 
(see Problem [T^) is strongly dual in the sense of combinatorial optimization. 
In particular, the vertex sets of the facets of P have a good characterization 
in terms of (in the sense of Edmonds [0). The problem is polynomial 
time equivalent to computing the cycles in Qp that correspond to the 2-faces 
of P. 

The problem can be solved in polynomial time in dimension at most three 
by computing a planar embedding of the graph, which can be done in linear 
time (Hopcroft and Tarjan Mehlhorn and Mutzel H). 

Related problems: |l|, |l^, |l| 

16. Facet System Verification for Simple Polytopes 

Input: The (abstract) graph Qp of a simple polytope P and a family T of 

subsets of nodes oi Qp 
Output: "Yes" if .F is the family of subsets of nodes oi Qp that correspond 

to the vertex sets of the facets of P, "No" otherwise 
Status (general): Open 
Status (fixed dim.): Open 

In j32| it is shown that both the "Yes"- as well as the "No" -answer can 
be proved in polynomial time in the size of Qp (provided that the integrity 
of the input data is guaranteed). Any polynomial time algorithm for the 
construction of an AOF- or US-orientation (see Problems |l^ and ^8|) of Qp 
would yield a polynomial time algorithm for this problem (see p^). 
Up to dimension three the problem can be solved in polynomial time (see the 
comments to Problems ^s] and jlj) . 

Related problems: |l|, |l^, |l|, |^ 
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17. AOF-Orientation 

Input: The (abstract) graph Qp of a simple polytope P 
Output: An AOF-orientation oi Qp 
Status (general): Open 
Status (fixed dim.): Open 

(Simple) polytopes admit AOF-orientations because every linear function in 
general position induces an AOF-orientation. 

In p2| it is shown that one can formulate the problem as a combinatorial 
optimization problem, for which a strongly dual problem in the sense of 
combinatorial optimization exists (see the comments to Problem hsh. Thus, 



the AOF-orientations of Qp have a good characterization (see Problem 16) in 
terms of Qp, i.e., there arc polynomial size proofs for both cases an orientation 
being an AOF-orientation or not (provided that the integrity of the input data 
is guaranteed). However, it is unknown if it is possible to check in polynomial 
time if a given orientation is an AOF-orientation. 

In dimensions one and two the problem is trivial. For a three-dimensional 
polytope P the problem can be solved in polynomial time, e.g., by producing 



a plane drawing of Qp with convex faces (see Tutte [62 ) and sorting the 
nodes with respect to a linear function (in general position). 
A polynomial time algorithm would lead to a polynomial algorithm for Prob- 
lem |6| (see [||). 

By duality of polytopes, the problem is equivalent to the problem of finding a 
shelling order of the facets of a simplicial polytope from the upper two layers 
of its face lattice. It is unknown whether it is possible to find in polynomial 
time a shelling order of the facets, even if the polytope is given by its en- 
tire face lattice. With this larger input, however, it is possible to check in 
polynomial time whether a given ordering of the facets is a shelling order. 

Related problems: |l^, |l^, ^ 



18. US- Orientation 

Input: The (abstract) graph Qp of a simple polytope P 
Output: A US-orientation oi Qp 
Status (general): Open 
Status (fixed dim.): Open 

Since AOF-orientations are US-orientations, it follows from the remarks on 
Problem |l^ that (simple) polytopes admit US-orientations and that the prob- 
lem can be solved in polynomial time up to dimension three. By slight adap- 
tions of the arguments given in [ p2[ , one can prove that a polynomial time 
algorithm for this problem would yield a polynomial time algorithm for Prob- 
lem ^ as well. 

In contrast to Problem no good characterization of US-orientations is 
known. 
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It is not hard to see that, by duality of polytopes, the problem is equivalent 
to the problem of finding from the upper two layers a partition of the face 
lattice of a simplicial polytope into intervals whose upper bounds are the 
facets (i.e., a partition in the sense of Stanley Similar to the situation 

with shelling orders, it is even unknown whether such a partition can be 
found in polynomial time if the polytope is specified by its entire face lattice. 
Again, with the entire face lattice as input it can be checked in polynomial 
time whether a family of subsets of the face lattice is a partition in that sense. 

Related problems: |l^, |l^, ^ 



4 Isomorphism 

Two polytopes Pi C TR"^^ and P2 C R'^^ are affinely equivalent if there is a 
one-to-one affine map T : aff(Pi) — > aff(P2) between the affine hulls of Pi 
and P2 with T{Pi) = P2- Two polytopes are combinatorially equivalent (or 
isomorphic) if their face lattices are isomorphic. It is not hard to see that 
afhne equivalence implies combinatorial equivalence. 

As soon as one starts to investigate structural properties of polytopes by 
means of computer programs, algorithms for deciding whether two polytopes 
are isomorphic become relevant. 

Some problems in this section are known to be hard in the sense that the 
graph isomorphism problem can polynomially be reduced to them. Although 
this problem is not known (and even not expected) to be AAP-complete, 
all attempts to find a polynomial time algorithm for it have failed so far. 
Actually, the same holds for a lot of problems that can be polynomially 
reduced to the graph isomorphism problem (see, e.g., Babai [^). 



19. Affine Equivalence of V-Polytopes 

Input: Two polytopes P and Q given in V-description 
Output: "Yes" if P is affinely equivalent to Q, "No" otherwise 
Status (general): Graph isomorphism hard 
Status (fixed dim.): Polynomial time 

The graph isomorphism problem can polynomially be reduced to the problem 
of checking the affine equivalence of two polytopes ||3^ . The problem remains 
graph isomorphism hard if 7i-descriptions are additionally provided as input 
data and/or if one restricts the input to simple or simplicial polytopes. 
For polytopes of bounded dimension the problem can be solved in polynomial 
time by mere enumeration of afhne bases among the vertex sets. 



Related problems: ^ 
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20. Combinatorial Equivalence of V-Polytopes 

Input: Two polytopes P and Q given in V-description 

Output: "Yes" if P is combinatorially equivalent to Q, "No" otherwise 

Status (general): coJ\fV-hard 

Status (fixed dim.): Polynomial time 

Swart j60) describes a reduction of the subset-sum problem to the negation 
of the problem. 

For polytopes of bounded dimension the problem can be solved in polynomial 
time (see Problems || and ^) . 

Related problems: ^, |l^, ^ 



21. PoLYTOPE Isomorphism 

Input: The face lattices Cp and Cq of two polytopes P and Q, respectively 

Output: "Yes" if £p is isomorphic to Cq, "No" otherwise 

Status (general): Open 

Status (fixed dim.): Polynomial time 

The problem can be solved in polynomial time in constant dimension (see 
Problem |2^ ). In general, the problem can easily be reduced to the graph 
isomorphism problem 



Related problems: 




22. Isomorphism of vertex- facet incidences 

Input: Vertex- facet incidence matrices Ap and Aq of polytopes P and Q, 
respectively 

Output: "Yes" if Ap can be transformed into Aq by row and column per- 
mutations, "No" otherwise 
Status (general): Graph isomorphism complete 
Status (fixed dim.): Polynomial time 

The problem remains graph isomorphism complete even if V- and 7i-descript- 
ions of P and Q are part of the input data . 

In constant dimension the problem can be solved in polynomial time by 
a reduction to the graph isomorphism problem for graphs of bounded 
degree, for which a polynomial time algorithm is known (Luks [^ ). 
Problem can polynomially be reduced to this problem. For polytopes of 
bounded dimension both problems are polynomial time equivalent. 



Related problems: ^ 
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23. Selfduality of Polytopes 

Input: Face Lattice of a polytope P 

Output: "Yes" if P is isomorphic to its dual, "No" otherwise 

Status (general): Open 

Status (fixed dim.): Polynomial time 

This is a special case of problem ^l[ In particular, it is solvable in polynomial 
time in bounded dimensions. 

It is easy to see that deciding whether a general 0/1-matrix A (not necessarily 
a vertex-facet incidence matrix of a polytope) can be transformed into A'^ 
by permuting its rows and columns is graph isomorphism complete. 

Related problems: |l| 
5 Optimization 

In this section, next to the original linear programming problem, we describe 
some of its relatives. In particular, combinatorial abstractions of the problem 
are important with respect to polytope theory (and, more general, discrete 
geometry) . We pick out the aspect of combinatorial cube programming here 
(and leave aside abstractions like general combinatorial linear programming, 
LP-type problems, and oriented matroid programming), since it has received 
considerable attention lately. 

24. Geometric Linear Programming 

Input: 7i-description of a polyhedron P C Q'', c G Q*^ 

Output: inf {c^x \ x G P} G QU{ — cx), oo} and, if the infimmn is finite, a 

point where the infimum is attained. 
Status (general): Polynomial time; no strongly polynomial time algorithm 
known 

Status (fixed dim.): Linear time in rri (the number of inequalities) 

The first polynomial time algorithm was a variant of the ellipsoid algorithm 
due to Khachiyan fssj] . Later, also interior point methods solving the problem 
in polynomial time were discovered (Karmarkar |^7| ). 

Megiddo found an algorithm solving the problem for a fixed number d of 
variables in 0{m) arithmetic operations (Megiddo [Q). 
No strongly polynomial time algorithm (performing a number of arithmetic 
operations that is bounded polynomially in d and the number of half-spaces 
rather than in the coding lengths of the input coordinates) is known. In 
particular, no polynomial time variant of the simplex algorithm is known. 
However, a randomized version of the simplex algorithm solves the prob- 
lem in (expected) subexponential time (Kalai Matousek, Sharir, and 
Welzl @). 

Related problems: ^ 
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25. Optimal Vertex 

Input: 7Y-description of a polyhedron P C Q'', c £ Q'' 

Output: inf {c^v \ v vertex of P} G Q U {00} and, if the infimum is finite, 

a vertex where the infimum is attained. 
Status (general): Strongly AAP-hard 
Status (fixed dim.): Polynomial time 

Proved to be strongly A/'P-hard by Fukuda, Liebling, and Margot By 
linear programming this problem can be solved in polynomial time if P is a 
polytope. In fixed dimension one can compute all vertices of P in polynomial 
time (see Problem]^). 

Related problems: |l|, ^ 



26. Vertex with specified objective value 

Input: 7i-description of a polyhedron P C Q'', c G Q"*, C G Q 
Output: "Yes" if there is a vertex v oi P with c^v — C; "No" otherwise 
Status (general): Strongly A/'P-complete 
Status (fixed dim.): Polynomial time 

Proved to be TV'P-complete by Chandrasekaran, Kabadi, and Murty and 
strongly A/'T^-complete by Fukuda, Liebling, and Margot The problem re- 
mains strongly A/'P-complete even if the input is restricted to polytopes |^ . 

Related problems: ^ 



27. AOF Cube Programming 

Input: An oracle for a function a : {0, 1}'' — > {+, — j"^ defining an AOF- 

orientation of the graph of the d-cube 
Output: The sink of the orientation 
Status (general): Open 
Status (fixed dim.): Constant time 

The problem can be solved in a subexponential number of oracle calls by 

the random facet variant of the simplex algorithm due to Kalai |36| . For a 

derivation of the explicit bound — 1 see Gartner . 

In fixed dimension the problem is trivial by mere enumeration. 

The problem generalizes linear programming problems whose sets of feasible 

solutions are combinatorially equivalent to cubes. 



Related problems: 
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28. USO Cube Programming 

Input: An oracle for a function a : {0,1}'' — > {+7~}'' defining a US- 
orientation of the graph of the d-cube 
Output: The sink of the orientation 
Status (general): Open 
Status (fixed dim.): Constant time 

Szabo and Welzl pi] describe a randomized algorithm solving the problem 
in an expected number of Oio!^^ oracle calls with a = •\/43/20 < 1.467 and a 
deterministic algorithm that needs 0(1.61'') oracle calls. Plugging an optimal 
algorithm for the three-dimensional case (found by Giinter Rote) into their 
algorithm, Szabo and Welzl even obtain an 0(1.438'') randomized algorithm. 
The problem not only generalizes Problem but also certain linear com- 
plementary problems and smallest enclosing ball problems. 
In fixed dimension the problem is trivial by mere enumeration. 

Related problems: ^ 
6 Realizability 

In this section problems are discussed which bridge the gap from combina- 
torial descriptions of polytopes to geometrical descriptions, i.e., it deals with 
questions of the following kind: given combinatorial data, does there exist 
a polytope which "realizes" this data? E.g., given a 0/1-matrix is this the 
matrix of vertex-facet incidences of a polytope? The problems of computing 
combinatorial from geometrical data is discussed in Section ||. 

The problems listed in this section are among the first ones asked in (mod- 
ern) polytope theory, going back to the work of Steinitz and Radermacher in 
the 1930's [|§. 

29. Steinitz Problem 

Input: Lattice L 

Output: "Yes" if L is isomorphic to the face lattice of a polytope, "No" 

otherwise 
Status (general): TVAT'-hard 
Status (fixed dim.): A/^T'-hard 

If L is isomorphic to the face lattice of a polytope, it is ranked, atomic, and 
coatomic. These properties can be tested in polynomial time in the size of L. 
Furthermore, in this case, the dimension d of a candidate polytope has to be 
rank£ — 1. 

The problem is trivial for dimension d < 2. Steinitz 's Theorem allows to solve 
c? = 3 in polynomial time: construct the (abstract) graph G, test if the facets 
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can consistently be embedded in the plane (linear time [ po|,p^ ) and check for 
3-connectedness (in linear time, see Hopcroft and TarjanpS]). 
Mnev proved that the Steinitz Problem for c?-polytopes with d + A vertices is 
A/'P-hard Even more, Richter-Gebert Q proved that for (fixed) d > 4 
the problem is A/'P-hard. 

For fixed d > 4 it is neither known whether the problem is in AfV nor whether 
it is in coAfV. It seems unlikely to be in MV, since there are 4-polytopes 
which cannot be realized with rational coordinates of coding length which is 
bounded by a polynomial in \C\ (see Richter-Gebert p3|). 

Related problems: ^ 



30. SiMPLiciAL Steinitz Problem 

Input: Lattice L 

Output: "Yes" if L is isomorphic to the face lattice of a simplicial polytope, 

"No" otherwise 
Status (general): TVAT'-hard 
Status (fixed dim.): Open 

As for Problem C is ranked, atomic, and coatomic if the answer is "Yes." 
In this case, the dimension d of any matched polytope is rank£ — 1. 
As for general polytopes (Problem |29|), this problem is polynomial time solv- 
able in dimension d < 3. 

The problem is TVP-hard, which follows from the above mentioned fact that 
the Steinitz problem for d-polytopes with d -I- 4 vertices is A/'P-hard and a 
construction (Bokowski and Sturmfels 0) which generalizes it to the sim- 
plicial case (but increases the dimension). It is, however, open whether the 
problem is A/"P-hard for fixed dimension. For fixed d > 4, it is neither known 
whether the problem is in MV nor whether it is in coAfV. 

The following question is interesting in connection with Problem |l^ (see also 
the notes there): Given an (abstract) graph G, is G the graph of a simple 
polytope? If we do not restrict the question to simple polytopes the problem 
is also interesting. 

Related problems: |l^, ^ 



7 Beyond Polytopes 

This section is concerned with problems on finite abstract simplicial com- 
plexes. Some of the problems listed are direct generalizations of problems on 
polytopes. Most of the basic notions relevant in our context can be looked up 
in for topological concepts like homology we refer to Munkres' book [flSj . 

A finite abstract simplicial complex Z\ is a non-empty set of subsets (the 
simplices or faces) of a finite set of vertices such that F G A and G C F imply 
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G G Z\. The dimension of a simplex _F e Z\ is \F\ — 1. The dimension dini(Z\) 
of A is the largest dimension of any of the simplices in Z\. If all its maximal 
simplices with respect to inclusion (i.e., its facets) have the same cardinality, 
then A is pure. A pure rf-dimensional finite abstract simplicial complex whose 
dual graph (defined on the facets, where two facets are adjacent if they share 
a common [d — l)-face) is connected is a pseudo-manifold if every (d — 1)- 
dimensional simplex is contained in at most two facets. The boundary of 
a simplicial {d + l)-dimensional polytope induces a c?-dimensional pseudo- 
manifold. 

Throughout this section a finite abstract simplicial complex A is given by 
its list of facets or by the complete list of all simplices. In the first case, the 
input size can be measured by n and m, the numbers of vertices and facets. 

31. EuLER Characteristic 

Input: Finite abstract simplicial complex A given by a list of facets 

Output: Euler characteristic x(/\) € TL 

Status (general): Open 

Status (fixed dim.): Polynomial time 

It is unknown whether the decision version "x(^) = 0?" of this problem is 
in MV . The problem is easy if A is given by a list of all of its simplices. For 
fixed dimension, one can enumerate all simplices of A and compute the Euler 
characteristic in polynomial time. 

Currently the fastest way to compute the Euler characteristic is to determine 
V = {5* : S' is an intersection of facets of Z\} and then compute x(^) in time 
0(|Vp) by a Mobius function approach, see Rota Q. Usually V is much 
smaller than the whole face lattice of Zi. V can be listed lexicographically by 
an algorithm of Ganter [ p3[ , in time ©(minjm, n} • a ■ |V|), where a is the 
number of vertex-facets incidences. 

Related problems: ^ 

32. /-Vector of Simplicial Complexes 

Input: Finite abstract simplicial complex A given by a list of facets 

Output: The /-vector of A 

Status (general): ^P-hard 

Status (fixed dim.): Polynomial time 

If A is given by all of its simplices the problem is trivial. Clearly, for fixed 
k, the first k entries of the /-vector can be computed in polynomial time, 
since the number of /c-simplices in A is polynomial in n. Hence the problem 
is polynomial time solvable for fixed dimension dim(Z\). 
It is unknown whether the decision problem "Given the list of facets of A 
and some (/? G N; is (^5 the total number of faces of Z\?" is contained in J\fV. 
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This problem is only known to be in MV for partitionable (see Problem p8| ) 
simplicial complexes (see Kleinschmidt and Onn ]39|]). 

To the best of our knowledge, no proof of #7^-hardness of the general problem 
has appeared in the literature. Therefore we include one here. 
Consider an instance of SAT, i.e., a formula in conjunctive normal form (CNF- 
formula) CiA- • -ACm with variables xi, . . . ,Xn (each Ci contains only disjunc- 
tions of literals). It is well known (Valiant |]63| ) that computing the number of 
satisfying truth assignments is ^P-complete. Define E — {ti, fi, . . . , tn, /„}. 
Part I. First, let E be the vertex set of a simplicial complex A defined by the 
minimal non- faces (circuits) C[, . . . , C^, Pi, . . . , P„, where Pi = {ti, fi} for 
every i. Here for any clause C, C :— {fj : Xj literal in C} U {tj : Xj literal 
in C}, e.g., for C ~ xi \/ X2 V X3 we have C = {/i, /2, ts}. The idea is that 
ti corresponds to the assignment of a true-value and fi corresponds to the 
assignment of a false- value to variable Xi. The circuits exclude subsets of E 
which include both ti and fi for all variables Xi and exclude truth-assignments 
to variables which would not satisfy a clause Cj. It is, however, allowed that 
for some variable Xi neither ti nor fi is included in a face. But every (n — 1)- 
face (n-subset of E) (if there exists any) corresponds to a truth-assignment 
to the variables (which uses exactly one value for each variable) and satisfies 
the formula. These subsets are counted by /„_i(Z\). Hence computing /„_i 
is ^P-complete and computing the /-vector of A is #P-hard. Moreover this 
shows that computing the dimension of a simplicial complex given by the 
minimal non- faces is AfV-haid. 

Part II. We now construct a simplicial complex A (the dual complex) which 
is given by facets. Define A by the facets C[, . . . , C^, Pi, . . . , P„, where for 
S E,S -.^ E\_S. We have that a set S Q_E is a face of A if and only if S 
is not a face of A. Hence, /„_i(Z\) -|-/„_i(Z\) — (^^). It follows that one can 
eflticiently compute /„_i(Z\) from /„_i(Z\). 

Related problems: |l| 



33. Homology 

Input: Finite abstract simplicial complex A given by a list of facets, i £ N 
Output: The i-th homology group of A, given by its rank and its torsion 

coefficients 
Status (general): Open 
Status (fixed dim.): Polynomial time 

There exists a polynomial time algorithm if A is given by the list of all 
simplices, since the Smith normal form of an integer matrix can be computed 
efficiently (Iliopoulos ^l|). For fixed i or dim(Z\)— i, the sizes of the boundary 
matrices are polynomial in the size of A and the Smith normal form can again 
be computed efficiently. 

Related problems: ^ 
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34. Shellability 

Input: Finite abstract pure simplicial complex A given by a list of facets 
Output: "Yes" if A is shellable, "No" otherwise 
Status (general): Open 
Status (fixed dim.): Open 

Given an ordering of the facets of A, it can be tested in polynomial time 

whether it is a shelling order. Hence, the problem in MV. 

The problem can be solved in polynomial time for one-dimensional complexes, 

i.e., for graphs: a graph is shellable if and only if it is connected. Even for 

dim(Z\) — 2, the status is open. In particular, it is unclear if the problem can 

be solved in polynomial time if A is given by a list of all simplices. 

For two-dimensional pseudo-manifolds the problem can be solved in linear 

time (Danarj and Klee [p^). 

Related problems: |l^, ^ 



35. Partitionability 

Input: Finite abstract simplicial complex A given by a list of facets 
Output: "Yes" if A is partionable, "No" otherwise 
Status (general): Open 
Status (fixed dim.): Open 

As in Problem p^, p artitionability is meant in the sense of Stanley |58| (see 



also [^). Noble |50 proved that the problem is in AfV. 
Partitionability can be solved in polynomial time for one-dimensional 
complexes, i.e., for graphs: a graph is partitionable if and only if at most one 
of its connected components is a tree. For two-dimensional complexes the 
complexity status is open. In particular, it is unclear if the problem can be 
solved in polynomial time if A is given by a list of all simplices. 

Related problems: |l|, || 
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